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ROCKET THRUST VARIATION WITH FOAMED LIQUID EROPEILANTS 

By G. Morrell 


SIMMARY 

Flow theory of liquid foams was applied to the rocket-engine injec- 
tion process. By foaming the propellants and thereby changing their hulk 
densities it is possible,, in theory, to vary rocket thrust continuously. 
An analysis of the method, assuming constant orifice flow coefficients, 
is presented and discussed. 

. Data from preliminary experiments in a 1000-pound-thrust ammonia - 

-f nitric acid rocket engine agreed only qualitatively with theory; two to 

=5 six times the theoretical gas -flow rate was required in the experiments. 

" > - It was demonstrated, however, that the "foam-flow" method of thrust vari- 

ation is feasible. 



INTRODUCTION 

Continuous thrust variation in rocket engines (throttling) is es- 
pecially desirable for piloted aircraft applications. In guided missiles 
the use of a single powerplant for both boost and sustainer operations 
would be possible; vernier thrust control near cut-off is another possi- 
ble field of application. Currently, rocket -thrust variation is accom- 
plished by using multiple fixed-thrust cylinders. 

Mechanical methods for varying thrust by injection orifice control 
are described in references 1 and 2. This report describes a method for 
controlling thrust by varying propellant density. Continuous density 
variation is achieved by foaming the propellants with an inert, insoluble 
gas. To maintain high propulsive efficiency, the exhaust -nozzle area 
must also be varied. This complementary problem of exhaust-nozzle-area 
control is treated in reference 3. 


Foam-flow theory is discussed in reference 4, and has been applied 
to the hydroduct, an underwater propulsion device, as reported in refer- 
ences 5 and 6 . Reference 6 also reports experimental data which are in 
substantial agreement with theory. 


UNCLASSIFIED 





2 


NACA RM E56K27 


In this report the foam-flow theory is applied to the injector ori- 
fice of a rocket engine to calculate the reduction in the propellant flow 
as a function of the ratio of the flow rates of the inert gas and the 
propellant. Other parameters involved are the initial injection pressure 
ratio, the gaa specific heat ratio, and a term including the fluid 
properties . 

Several preliminary experiments were conducted in a nominal 1000- 
pound-thrust ammonia - nitric acid rocket to measure the reduction in 
liquid flow over a small range of gas flows. The deviation of the re- 
sults from calculated values is discussed as well as the feasibility of 
this method of thrust control. 


SYMBOLS 

The following symbols are used in this report: 

A cross-sectional area 

Cp orifice discharge coefficient 

g conversion constant 

p pressure 

R gas constant 

T absolute temperature 

u velocity 

V volume 

W weight 

w weight-flow rate 

a ratio of injection pressure to chamber pressure for full flow, 
p^/p c q , dimensionless (see subscript list) 

0 RTp^/p c,o> <3iniens ionless 

X ratio of specific heats, dimensionless ■ 

ratio of part-thrust flow rate to f ull -thrust flow rate, w^/w^ Q , 
dimens ionles s 
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t ratio of gas-flow rate to liquid-flow rate, Wg/w^, dimensionless 
p propellant density 

t ratio of part -thrust burning time to total burning time, 
dimensionles s 

CD ratio of discharge coefficients, Cj^ Cj),f, dimensionless 
Subscripts: 

c combustion chamber 

f foam 

g gas 

I liquid 

0 full -thrust condition 

th theoretical 

1 upstream of injection orifice 

2 injection orifice discharge station 


ANALYSIS 

In the following analysis, expressions are derived for the reduction 
in liquid flow as a function of the gas 'to liquid ratio s for several 
values of a and p and for "A Vhlue of y of 1.67. Two flow eases are 
treated: compressible, isothermal gas flow, assuming thermal equilibrium 

between the gas and the liquid and essentially constant internal energy; 
and compressible adiabatic flow, assuming no interchange of energy between 
the gas and the liquid. In both cases, any heat transfer between the 
fluid and the surroundings is ignored. 


Derivation of Flow Equations 

The foam-flow equations are derived on the basis that the gas is 
uniformly dispersed and is insoluble in the liquid and that the flow is 
one-dimensional. The continuity equation is 

w g + w 2 = W 2 (l + e) = pAu 


CD 
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Gas and liquid volumes may "be considered additive, so that 

V - Vg + V Z 

or since 

v w 
V a - 

P 


w g + n w g . wi 

m -■■■ ■ *r — 

p p g p i 

Substituting the gas -liquid ratio 


s 




gives the following equation for foam density: 


P 


e + 1 



( 2 ) 


The simplified Euler equation for one -dimensional compress ible-f low is 
written as 


+ &. 0 (3) 

g p s 

In solving equation (3), the following assumptions were made: (l) flow 

velocity at station 1 Is negligible compared with that at station 2; 

(2) at station 2 the pressure is equal to the combustion-chamber static 
pressure, that is, p 2 = p c ; (3) the characteristic velocity based on 

liquid flow (p c A^.g/w^, where Aj. is nozzle-throat area) remains constant, 
or p c oc w-j; (4) the liquid flow is isothermal; (5) the gas behaves ideal- 
ly, that is, Pg “ p/RT; and (6) the propellant mixture ratio remains 
constant . 

Integration of equation (3) for the case of isothermal gas flow 
after substituting equation ( 2 ) yields 



Pi - ?c \ 

pi / 


^2 


(4) 
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and substitution in equation (l) gives 


_ 90,^2 J 2g / 

A 2 “ l + «fl+«\ 


„ Pi Pi - P c ' 

6R!T In — + 

P c P l j 


c D,f Yttt ^ ^ + 

SET 1 


Pi Pi - Pc ' 

Pc Pi / 


Pc p 7 


where the discharge coefficient has been introduced. For no gas flow, 
6 = 0, wj — v l#Q , P c =.p C;0 , and 

“j^T “ C D,I V 2 S Pl(Pl " P c ,o) , 0 

Dividing equation (5) by (6) gives the liquid-flow ratio. 


1 P 
+ 6 L 


6RTp x In Pl /p c + (p-l - p c )‘ 


Pi - Pc.o 


^sRTP^ 


By letting p, /p = a, and KCp,/p = 8, and since 

1* C j O * C;0 * 

5 = = P c /p c , 0 > equation (7) may be reduced to 


^ m | + SL 

1 + e yx - 1 8 a 

“(if + A 


e 2 p 2 +• 2ep6 + 8 2 


ep In ^ + a - 5 
ai 2 (l +e) (a - 1) 
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If on the other hand the gas -phase flow through the Injector is 
adiabatic, the gas density derived from the ideal gas law and the adia- 
batic equation of state is 


p^Ap-]_ Ct~i )/t 

p g " 

which upon substitution with equation (2) in equation (3) followed by 
integration yields 



( 9 ) 


Following the same procedure as above and substituting equation (9) 
in equation (l) gives the liquid flow per unit area 



Division of equation (10 ) by w^o/Ag 
liquid ratio 


from equation (6) yields the 





r r=i 

r=y 




sr RIIl i 

Pl r 

. p r 

■^c 

Pi - Pc 

v l 1 

S3 — 

/ 1+ ‘ 

^ H 
Pi 

1 


1 

o 

o 

1 

p l< p l - p <=,o> 


v l 


,o 



( 11 ) 
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Again, -by letting 5 s - vJv C)0 , a » Pi/p Cj0 , and 

P =* HEjPj/p Q , equation (ll) becomes 




1 

/ni ni\ 

eyP / a T - 5 7 

ai 2 (l+s) 

r-l\ 0,-1 

LCr-iM V / 



( 12 ) 


Equations (8) and (12) were solved on a card-programmed calculator 
for several values of a, and p, for a single value of x (helium gas 

assumed; x> 1*67), and for a constant discJiarge coefficient, that is, 

when <x> equals 1. The last condition was assumed in order to simplify 
the calculation. There is no experimental evidence either to support or 
negate this assumption. Figures 1 and 2 show the relations between the 
liquid ratio 5 and the gas-to-liquid ratio e for compressible, iso- 
thermal flow and compressible, adiabatic flow, respectively. The depend- 
ence of s on a is shown in figure 3, and the dependence of s on p 

is shown in figure 4 both for a 5 value of 0.5. 


Discussion 

In terms of gas requirements, adiabatic gas flow is less economical 
than isothermal flow. From the nature of the flow process, especially 
the short stay time in the injector, it appears probable that the amount 
of energy interchange between liquid and gas will be small, and the gas- 
phase flow will be more nearly adiabatic than isothermal. 

For the a range normally encountered in rocket practice (1.1 to 
1.3), injection pressure has a large effect on gas economy. From figure 
3 it appears that engines using the foam-flow technique should be de- 
signed with injection pressure ratios in the range of 1.5 to 2.0, if 
possible. 
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The parameter 8 also has a large influence on gas economy; the 
larger the value of the lower is the value of 6 for a given value 
of 8, as shown in figure 4. In effect, this means that for a given gas 
and chamber pressure, foam-flow thrust variation is most practical for 
high-density propellants. In the following table are listed values of 
8 for several propellants for the case where R is 386 foot pounds per 
pound per °R (helium gas) and p c 0 is 450 pounds per square inch; 


Propellant 

* 1 > 

°R 

Pi 

lb/cu ft 

P 

Hydrogen 

37 

4.4 

0,97 

Oxygen 

167 

71 

70.7 

Ammonia 

460 

41.3 

113.0 

Jet fuel 

530 

47 

150.7 

Nitric acid 

530 

93 

290.1 


It is apparent immediately that fuel - nitric acid systems are best 
suited to the foam-flow method of thrust variation, and that the method 
would be impractical for hydrogen systems. 

To illustrate the application of the analysis, the following example 
is presented. The model chosen is a rocket missile using jet fuel and 
nitric acid at a mixture ratio of 4. Full -thrust conditions were assumed 
to be as follows: thrust, 20,000 pounds; net specific impulse, 220 pounds 

per second per pound; chamber pressure, 500 pounds per square inch; in- 
jection pressure, 650 pounds per square inch; and injection temperature, 
530° R. Powered flight duration was held constant at 2 minutes, and 
helium gas was used, which was stored at 3000 pounds per square inch and 
530° R. For maximum tot el impulse the volume of propellants required is 
140 cubic feet. Helium density at storage conditions, 1.93 pounds per 
cubic foot, was calculated from the data of reference 7. The average 
density of the propellants is 78 pounds per cubic foot, equivalent to a 
8 value of 222. When compressible, adiabatic flow is assumed, the rela- 
tion of 8 and s from figure 2 is 


8l 

"0775 1 

0.5 1 

0.25 


0.0085 

0.0235 

0.065 


These data and assumptions were used to calculate the required gas and 
liquid volumes as a function of the fraction of total powered flight 
time for one-quarter, one -half, and three-quarters thrust operation. 

The results of the calculations are shown In figure 5. To provide half 
thrust for half the flight duration would require a gas volume equal to 
about 30 percent of the liquid volume. Total volume would be only 
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slightly less than that for full -thrust,, full-duration flight. Since 
high-pressure gas tankage will he heavier than liquid tankage on a unit 
volume basis, it is clear that the foam-flow technique for thrust varia- 
tion imposes a weight penalty on the vehicle structure, even though the 
gross vehicle weight is less than that for full-thrust flight. 

How this weight penalty compares with the added weight due to mul- 
tiple thrust cylinders and the gear for varying injection orifice area 
will depend on the flight program and vehicle design and cannot be stated 
in generalized terms. Qualitatively, however, it is apparent from figure 
5 that the foam-flow technique is most advantageous for smaller values of 
t and larger values of 5. For large values of 5, precise control of 
gas flow may be a problem. 


EXPERIMENTAL APPARATUS 

A short series of experiments were run in a nominal 1000 -pound- 
thrust water-cooled rocket to check the analysis. The propellant system 
was ammonia - white fuming nitric acid (WFNA) with helium used as the 
pressurizing and foaming gas. Each propellant was foamed separately 
before entering the injection manifold. Flow-line addition of lithium 
metal to the ammonia caused spontaneous ignition in the combustor. A 
schematic drawing of the flow system is shown in figure 6. 

The thrust cylinder had an inside diameter of 4 inches and an over- 

3 

all length of 14=- inches. The cylindrical portion of the combustor was 

O 

5 13 

8^ inches long; nozzle-throat diameter, 1^ inches; and nozzle-exit dism- 

3 

eter, 3^ inches. The characteristic length of the rocket was 47 inches. 

The doublet -type injector consisted of 24 pairs of fuel and oxidant 
orifices arranged in a circle with a mean diameter of 2.5 inches. The 
fuel-orifice diameter was 0.0515 inch, and the oxidant-orifice diameter, 
0.072 inch. 


Gas -Injection Device 

Several gas -injection devices were qualitatively tested with water 
in a transparent tube to check the degree of gas dispersion and flow 
stability. The final design, shown in the insert of figure 6, consisted 
of a 2-inch length of 5/l6-diameter tubing with 220 holes, each with a 
0.0135-inch diameter, arranged in 11 circles of 20 holes each. The tube 
was fitted into a T so that the tube axis coincided with the thorough axis 
of the T. 
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Calibration of the gas-injection device with water as the liquid In 
the propulsion system indicated that flow was free from surges when the 
gas-injection pressure was not more than 100 pounds per square inch 
greater than the liquid-injection pressure. 


Instrumentation 

Flow rates of the propellants were measured by rotating- vane-type 
meters with an error of about 1 percent. Gas -flow rates were measured 
by orifices with an error of about 2 percent. Pressures were measured 
with strain-gage -type transducers having an error of abgut 2 percent. 
Thrust was measured with a strain-gage load cell having an accuracy of 
about 1 percent. Temperatures were measured with chromel-alumel 
thermocouples . 


Operating Procedure 

The engine was started and operated at full-flow conditions for 
about 10 to 15 seconds, then helium was admitted and the engine operated 
at part thrust for an additional 10 to 15 seconds. Operation was finally 
shifted to full-flow conditions for approximately 3 to 5 seconds prior to 
shut down in order to check the initial data. 


RESULTS AND DISCUSSION 

Data for the four experiments performed are listed in table I. 

Table II shows the values of the analytic parameters and a comparison of 
the experimental and theoretical gas-liquid ratios. The assumption of 
constant characteristic velocity (chamber pressure proportional to flow) 
was nearly fulfilled in the experiments as shown in table I. The largest 
deviation occurred in the second run. The assumption of a constant mix- 
ture ratio was not fulfilled in the first two runs, (i.e., 5 for the fuel 
and the oxidant were not the same (table H)), which may account for some 
of the variation in the characteristic velocity. 


Comparison of Gas-Liquid Ratios 

Although the experiments definitely show that substantial thrust 
variation can be obtained by foaming the propellants to alter their bulk 
density, the experimental gas-liquid ratios were about twice the theoret- 
ical values for the fuel and three to six times the theoretical values 
for the oxidant. 

There are a number of probable reasons for the discrepancy between 
the experimental and predicted results: (l) the analysis assumed the 
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same discharge coefficient for foam and liquid, and the experiments de- 
* parted seriously from this condition, as is shown in table III (discharge 

coefficients were calculated by eqs. (6) and (10 ) ) j (2) the gas-injection 
device may not have produced, a uniform gas dispersion since the design 
was selected on the basis of -gross visual observations y and (3) the in- 
jector manifold design tended to produce a centrifugal component of ve- 
locity in the fluid which may have caused gas-liquid separations. The 
anomaly of discharge coefficients greater than unity may have been caused 
by experimental error. 

3 

j More research is needed to improve the foam-flow method of rocket- 

' thrust variation. Studies are needed on the methods for producing uni- 

form, stable foams and for the design of flow passages which will minimize 
gas-liquid separation. The flow of foams through orifices should also be 
studied to establish discharge coefficients as a function of the flow 
conditions and the fluid properties. Since one of the propellants will 
be used as a coolant, studies should be conducted on the heat-transfer 
characteristics of foamed liquids. It might be expected that for moder- 
ate gas -liquid ratios, the heat-transfer characteristics of foams will 
approach those for liquids under nucleate boiling conditions. It was 
found in water calibrations of the flow systems that large pressure dif- 
ferences between liquid and gas produced intermittent flow of liquid. 

Flow of this type would certainly result in low-frequency combustion os- 
cillations in a rocket engine. When applying the foam-flow technique to 
a powerplant, therefore, special attention will be required in the design 
of the flow system to avoid this type of flow. 


Operating Characteristics 

The transition from full-thrust to part-thrust and back to full- 
thrust was free from surge in all runs. Since the pressure difference 
between liquid and gas was maintained at a low value, the combustion was 
also free from oscillations. 


CONCLUDING REMARKS 

This work has demonstrated that it is technically feasible to vary 
rocket -engine thrust by the foam-flow method. Research on several prob- 
lems is needed in order to improve the efficiency of the process, partic- 
ularly in view of the fact that two to six times the theoretical gas-flow 
rates were required in the few experiments performed. Further than this, 
detailed analyses are required to establish the relative merits of foam- 
flow thrust variation against mechanical techniques such as orifice-area 
variation or pump-speed variation. It Is quite probable that no one 
method will be best, but rather that the application and propellants 
chosen will deter min e the optimum method. 
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SUMMARY OF RESULTS 

An analysis is presented, on a method for continuously varying rocket 
thrust by foaming the propellants to change their bulk density. For a 
given thrust ratio, the gas-liquid ratio depends on the full -thrust- 
injector pressure ratio, the molecular veight and temperature of the gas, 
the liquid density, the normal combustion pressure, and the gas expansion 
process. 

From a series of four experiments in a nominal 1000 -pound-thrust 
rocket using ammonia - nitric acid propellants -with helium as the foaming 
gaB, the following results were obtained: 

1. The technical feasibility of the foam-flow method for thrust 
variation was shown. 

2. The experimental gas -liquid ratios were two to six times the 
predicted values. 

3. The characteristic velocity for part-thrust operation was within 
10 percent of that for full-thrust operation, indicating that foaming of 
the propellants does not impair combustion efficiency. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 27, 1956 
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ON ROCKET THRUST VARIATION 


Oxidant 

Oxidant- 

Fuel- 

Oxidant- 

Fuel- 

Fuel- 

Oxidant- 

temper- 

helium 

helium 

in J ection 

injection 

helium 

helium 

ature, 

°F 

temper- 

ature, 

°F 

temper- 

ature, 

°F 

pressure, 
Ib/sq in* 

pressure, 
Ib/aq in. 

pressure, 
Ib/sq in. 

pressure, 
ih/sq in. 

86 

— 

— 

350 

550 

— - 

— 

86 

45 

48 

350 

560 

550 

350 

52 

— 

— 

350 

550 

— 

— 

52 

36 

42 

350 

550 

550 

350 

32 

— 

— 

350 

550 

— 

— 

52 

30 

36 

350 

550 

550 

400 

52 

— 

-- 

400 

550 

— 

— 

52 

28 

32 

400 

550 

550 

450 
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TABLE H. - COMPARISON OF EXPERIMENTAL AND 


THEORETICAL GAS -LIQUID RATIOS 


Run 

8 

Pc/Pc,o 

e 

— 

CL 

P 

Theoretical 
gas -liquid 
ratio for 
compress- 
ible, adia- 
batic flow 

A 

-P 

<0 

Oxidant system 

1 

0.720 

0.598 

0.00983 

1.45 

518 

0.0033 

2.98 

2 

.830 

.625 

.0122 

1.39 

500 

.0019 

6.42 

3 

.555 

.537 

.0344 

1.56 

515 

.0057 

6.04 

4 

.516 

.534 

.0367 

1.66 

484 

.0065 

5.65 

Fuel system 

1 

0.550 

0.598 

0.0174 

2.28 

226 

0.0097 

1.79 

2 

.494 

.625 

.0241 

2.19 

211 

.0131 

1.84 

3 

.565 

.537 

.0161 

2.29 

219 

.0094 

1.72 

4 

.595 

.534 

.0170 

2.16 

210 

.0092 

1.85 


TABLE III. - COMPARISON OF INJECTOR 
DISCHARGE COEFFICIENTS 


Run 

Full -thrust 
coefficient, 

^,2 

Part -thrust 
coefficient, 

°D,f 

Coefficient 

ratio, 

CD 

Oxidant system 

1 

0.507 

0.910 

0.558 

2 

.489 

1.02 

.479 

3 

.451 

1.08 

.418 

4 

.479 

1.1 

.435 

Fuel system 

1 

0.590 

0.886 

0.666 

2 

.611 

.875 

.699 

3 

.637 

.975 

.653 

4 

.617 

.968 

.638 | 
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(a) Ratio of injection pressure to chamber pressure for full flow, 1.1. 
Figure 1. - Variation of liquid ratio for isothermal gas flow. 
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CD 

m 

CVJ 



(b) Ratio o f injection pressure to chamber pressure 
for full flow j 1.3. 


Figure 1. - Continued. Variation of liquid ratio for isothermal gas flow. 
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0 * 01 .02 . 03 . 04 . 05 . 06 

Gas -liquid ratio, e 

(c) Ratio of injection pressure to chamber pressure 
for full flov, 1.9. 

Figure 1. - Concluded. Variation of liquid ratio for isothermal gas flow. 
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(a) Ratio of injection pressure to chamber pressure 
for full flow, 1.1. 


Figure 2. - Variation of liquid ratio for adiabatic gas flow with, a specific heat 
ratio of 1.67. 













Gas -liquid ratio 
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Initial Injection presBure ratio, a 

Figure 3. - Variation of gas-liquid ratio with initial Injection 
pressure, ratio; liquid ratio of 0.5. 




Parameter 
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Gas-liquid volume ratio Total volume of gas plus liquid, cu ft 
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0 .2 .4 .6 .8 

Ratio of part-thrust burning time to total 
burning time , t 


Figure 5. - Gas and liquid requirements for part- 
thrust flight. 
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Helium gas 



Foaming device 


Figure 6. - Schematic diagram of experimental rochet system. 
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